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Abstract We investigate the dynamics of two- component Bose-Josephson junction 
composed of atom- molecule BECs. Within the semiclassical approximation, the multi- 
degree of freedom of this system permits chaotic dynamics, which does not occur in 
single-component Bose-Josephson junctions. By investigating the level statistics of the 
energy spectra using the exact diagonalization method, we evaluate whether the dy- 
namics of the system is periodic or non-periodic within the semiclassical approximation. 
Additionally, we compare the semiclassical and full-quantum dynamics. 



1 Introduction 

Since its experimental realization, Bose-Einstein condensation in dilute atomic gases 
has been offering opportunities to research macroscopic quantum phenomena. Espe- 
cially, one of the most fascinating macroscopic quantum phenomena is the Joseph- 
son effect between two Bose-Einstein condensations(BECs) trapped in a double-well 
potential, which is called the Bose-Josephson Junction(BJJ). BJJ has been realized 
experimentally, which allowed direct observation of the dynamics of macroscopic wave- 
functions directly—. Though in BJJ the spatial coherence of BECs is focussed, the 
Josephson effect occurs not only between spatial separated BECs but also between 
internal degrees of freedom in a single BEC. In particular, the Josephson- like effects 
between atomic and molecular states have been discussed theoretically^^. In this pa- 
per, we consider atom-molecule BECs in a double- well potential. In this two component 
BJJ, the internal atom-molecule tunneling has significant influence on the ground state 
property and dynamics. As shown later, this internal tunneling induces the localized 
ground states. As for the dynamics, the multi-degree of freedom of this system permits 
chaotic dynamics, by adding the atom-molecule tunneling as a new degree of freedom. 

First, we explain the semiclassical dynamics of BJJ. In a single component BJJ, 
the dynamical degrees of freedom are the relative particle number and the relative 
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phase, so the dimension of the system is two. Since the trajectories in the phase space 
cannot intersect with each other, significantly complex trajectories cannot exist in 
2D systems. However, in higher dimensional systems, the trajectories in the phase 
space travel around much more freely than 2D dynamical systems. This leads to richer 
dynamical behaviors in a two component BJJ, and the chaotic dynamics can occur. 

In general full-quantum systems, although the chaotic dynamics can not occur(see 
Ref . 5 and references there in) , the quantum signature of the semiclassical non-periodicity 
appears in the level spacing distribution—. In our system, we can control the multiplic- 
ity of degrees of freedom by varying the atom-molecule internal tunneling strength. 
Based on these considerations, we investigate the statistical property of energy spectra 
and compare the semiclassical and full-quantum dynamics. 



2 Model and approximations 

In this section, we explain our model and approximations. The second- quantized Hamil- 
tonian for Bose atoms and molecules can be written as 

fc2 



i=a,b i=a,b 

+g ab J <h$$t$b$ a ~ X J dr (*b*a*a + # 6 #M) + <5 j dr^ b W b , (1) 

where ifr a and represent field operators for Bose atoms and molecules respectively, 
A is the coupling strength between atomic and molecular states, 5 is the energy differ- 
ence between atoms and molecules, and V e ^ct(r) is the double- well potential. The inter- 
atomic, the inter-molecule, and the atom-molecule interactions can be approximated 
in terms of the s-wave scattering lengths as gi — 4.-K?i 2 a s i/mi,g a b — 67r^ 2 a sa 5/m a (« = 
a^b^rrib — 2m a )- Here, m a is the mass of a Bose atom. Furthermore, we introduce the 
four-mode approximation, which concentrates on condensate modes only, and neglect 
the effect of the particles occupying other modes. Then, field operators can be approx- 
imated as # a ~ $aLO>L+&aRa>R,&b - ®bL b L + ®bR b R, where $ a L,®aR(®bL,®bR) are 
the wavefunctions of the atomic (molecular) condensate modes in the left well and the 
right well respectively, ac, 6r) are annihilation operators for the atomic (molecular) 

condensate modes in the left well and the right well respectively. Applying these ap- 
proximations to Eq.(l), we obtain the quantum four-mode Hamiltonian (four- mode 
model). 



H = -J a (a[a R + <4a L ) - J b (b[b R + b ] R b L ^ + A (b ] L b L + b^b^ 

+ lf ( a L a L a ^ a ^ + a R a R a R a R) + Y { b L b L b L b L + b R b R b R b R) 
+U ah {a ] L a L b ] L b L + a ] R a R b ] R bR^j 

-g (b^ L a L a L + b^ R a R a R + b L a\a\ + b R a^ R a^ . (2) 
Here, the parameters are defined as follows : 
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Ui = g 

U a b = dab 



h 2 9 

-^V 2 + Kxt(r) 
Zrrij 
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2m, 



V z + T4xt(r) 



,L\ 2 \$ bL \ 2 , 



J dr\<P iL \ 4 = 9i J dr\<P iR \ 4 , 

J dr\$ aR \ 2 \$ bR \ 2 = g ab J dv\$ al 
g = X J dr$Z L $ aL $ aL = \ J dv$* bR $ aR $ aR , 
A = 5 f dv\$ bL \ 2 + El -2E° a = sJ dr\$ bR \ 2 + E° b - 2E° a , 



$iR, (4) 
(5) 
(6) 
(7) 
(8) 



where % — a represents atomic BEC modes, i = b represents molecular BEC modes, 
and L,R express the left well and right well respectively. And also, we use the notation 
as Za — (o) l cll - a R dRj /N as needed. 



3 Results 

First, we investigate the effect of the atom-molecule tunneling on the ground state. 
Next, we will show the relation between the level statistics and the periodicity of dy- 
namics. We choose the parameters as follows. The ratio A = NU a /(2J a ) is estimated 
as 15 in the single-component experiment 1 , and the atomic interaction strength nor- 
malized by the atomic tunneling strength can be obtained as U a /Ja — 3 x 10 -2 . We 
use this value for the atomic interaction strength. As for the molecular interaction 
strength, we suppose that the molecular scattering length is the same as the atomic 
one and that the shape of condensate wavefunctions of atoms and molecules are also 
the same. Under this condition — U a /2 from Eq. J5}. In this study we set the total 
particle number as TV = 20. We next consider the atom-molecule interaction. The ex- 
periment 7 indicates that the atom-molecule interaction of 87 Rb is attractive, and thus 
we choose as U a b/J a — —2.3 x 10 -2 . Finally, we set A/J a = — 1. 



3.1 Phase transition of the ground state 

Within the semiclassical approximation, the particle localization transition in the 
ground state is induced by the atom- molecule internal tunneling-. The ground state 
at y/Ng/Ja < 2.55 is symmetric, which has equal particle populations in each well. At 
y/Ng/ Ja — 2.55, the ground state become degenerate. For y/~Ng/J a > 2.55, the ground 
state obtains non-equal particle population. Here, we investigate this transition within 
a full-quantum treatment. By the exact diagonalization of the four- mode Hamiltonian 
J2]), we obtain the eigenstates and energy eigenvalues. 

The ^-dependence of (\z a \) is presented in Fig.[T] (\z a \) increases rapidly at V~Ng/ J a — 
3.5. From Fig. [2] we see that the fluctuations of \z a \ is maximum at this point. Further- 
more, the ground state become degenerate for y/~Ng/J a > 3.5, as shown below. The 
^-dependence of excitation energies is shown in Fig. [3] Increasing the atom-molecule 
tunneling g, the first excitation energy goes to zero at y/Ng/Ja — 3.5. For larger g, the 
ground state and some low lying energy levels become degenerate. These facts indicate 
that a phase transition to the localized ground state occurs at y/Ng/J a — 3.5. Finally, 
we comment on the system-size dependence of the results. We find that the results for 
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Fig. 2 (color online) (A\z a \) of N = 20(red 
Fig. 1 (color online)The ^-dependence of points) and TV = 30(green points). (Z\|z a |) is 

denned as (A\z a \) = y/(zl) -(\z a \) 2 . 




Fig. 3 (color online)g-dependence of the energy spectra normalized by J a - Lower 30 levels 
are presented. 



(\z a \) are qualitatively the same for different values of N. The size dependence appears 
more prominently in fluctuations; as shown in Fig. [2] the peaks of (Z\|z a |) at the critical 
point are sharpened for larger N. This indicates the occurrence of the quantum phase 
transition in the thermodynamic limit. 



3.2 Level statistics and dynamics 

Within a semiclassical approximation, we found the existence of non-periodic dynam- 
ics—. This motivates us to investigate the full-quantum counterpart. 

In this section, we investigate the effect of the atom-molecule tunneling on the 
periodicity of dynamics. The semiclassical dynamics is periodic in the small and large 
VNg/Ja region, however, non-periodic in the intermediate region. Considering this by 
full-quantum treatment, we show how the level spacing distributions change by varying 
the strength of the atom-molecule tunneling. From level spacing distribution, one can 
decide whether the dynamics of the system is periodic or non-periodic 6 . 
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Fig. 4 (color online)Level spacing distribu- 
tions at \fNg/J a = 0.447. Blue and violet 
bars represent odd and even parity energy- 
spectra. Red and green lines are Poisson and 
Wigner distributions. 
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r = (J a /h)t 

Fig. 5 (color online)The time evolution of 
N aL = (a ] L a L ) at VNg/J a = 0.447. As the 
initial condition, we shift the atomic particle 
distribution of the ground state to the left well 
by a single atom. Green and red lines repre- 
sent semiclassical and full-quantum time evo- 
lution. 



When a system exhibits regular motions in a semiclassical approximation, level 
spacing distribution shows Poisson-type distribution P(S) = e~ S except some cases 10 , 
where S is the level spacing normalized by mean level spacing. On the other hand, when 
the semiclassical dynamics of the system is chaotic, level spacing distribution is Wigner- 
type distribution P(S) — (vr/2) S exp (— 7rS 2 /4). In this section, we compare the level 
spacing distribution of the four-mode Hamiltonian J2) with Wigner and Poisson dis- 
tribution in order to investigate the periodicity of semiclassical dynamics. In addition, 
we show that the time evolutions in a semiclassical approximation and a full-quantum 
treatment 11 are different in a chaotic regime. Before the analysis, we need to classify 
the energy spectra according to the symmetry of the system—. Since the double-well 
potential has the left-right symmetry, we should investigate the odd and even parity 
spectra separately. 

From Fig. [4] the level spacing distiribution in the small g regime is very close 
to Poisson distribution. In this regime, Fig. \E\ shows that the time evolutions in a 
semiclassical approximation resemble the full-quantum dynamics, and both exhibit 
regular motions. Furthermore, the level spacing distribution is Poisson-like in the larger 
g region(see Fig. [6}. The time evolutions of N a L = (a^a^) exhibit regular motion in 
both semiclassical and full-quantum treatments, as shown in Fig. [7] In contrast to Fig. 
El we see a remarkable difference between the semiclassical and full-quantum dynamics. 
This is due to the particle localization in one well, which decreases the particle number 
in another well and enhances quantum fluctuations. In VNg/Ja ^> 1 and VNg/Ja <C 1 
regions, internal or interwell tunneling are dominant, and thus the multi-degree of 
freedom of this system does not have strong influence. In contrast, as shown in Fig. [8) 
the almost all energy spectra obey Wigner distribution in the intermediate coupling 
region. Fig. [9] shows that the semiclassical dynamics is chaotic and quite different from 
the full-quantum time evolution. It is shown that the quantum fluctuations induce 
dynamical localization. The similar behavior was found in a kicked roto r 13 14 . 

4 Conclusions 

In this paper, we showed the particle localization in ground states as shown in Fig. 
[TJ Furthermore, we investigated how the internal tunneling deform the level spacing 
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Fig. 6 (color online)Level spacing distribu- 
tions at \f~Ngl J a = 33.5. Blue and violet bars 
represent odd and even parity energy spectra. 
Red and green lines are Poisson and Wigner 
distributions. 
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r = (J a /h) t 

Fig. 7 (color online)The time evolution of 
N aL = (a ] L a L ) at y/Ng/J a = 33.5. As the 
initial condition, we shift the atomic particle 
distribution of the ground state to the right 
well by a single atom. Green and red lines 
represent semiclassical and full-quantum time 
evolution. 
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Fig. 8 (color online)Level spacing distribu- 
tions at y/Ng/Ja = 3.43. Blue and violet bars 
represent odd and even parity energy spectra. 
Red and green lines are Poisson and Wigner 
distributions. Lower 10 % and higher 10 % 
levels are omitted. 
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r = (J a /h) t 

Fig. 9 (color online)The time evolution of 
N aL = (a ] L a L ) at VNg/J a = 3.43. As the 
initial condition, we shift the atomic particle 
distribution of the ground state to the right 
well by five atoms. Green and red lines repre- 
sent semiclassical and full-quantum time evo- 
lution. 



distribution and dynamics. Increasing the internal tunneling, the dynamics change 
to be non-periodic. Further increasing the internal tunneling, periodic time evolution 
emerges again. In addition, we showed that semiclassical and full-quantum dynamics 
is quite different in chaotic region; dynamical localization occurs within full-quantum 
treatment. 
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